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Abstract 

We apply the method of unitary transformations to a model two nucleon potential and 
construct from it an effective potential in a subspace of momenta below a given cut-off A. 
The S-matrices in the full space and in the subspace are shown to be identical. We solve 
numerically the Schrodinger equation in the small momentum space and recover exactly the 
bound and scattering states of the full theory. We then expand the heavy repulsive meson 
exchange of the effective potential in a series of local contact terms and discuss the question 
of naturalness of the corresponding coupling constants. Using our exact effective theory we 
address further issues related to the chiral perturbation theory approach of the two-nucleon 
system. The coordinate space representation of the effective potential is also considered. 
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1 Introduction 



Chiral perturbation theory is hoped to provide insight into the nuclear force problem and pos- 
sibly even lead to a quantitative framework. Chiral symmetry imposes constraints on possible 
momentum and spin dependences of the nuclear forces but the framework is restricted to mo- 
menta below a certain scale. Only in this regime one can set up a power counting scheme which 
limits the number of interaction terms in the nuclear Hamiltonian. Whether this scheme works 
quantitatively and is applicable to interacting nucleons seems not yet fully established, despite a 
large number of investigations and considerable progress over the last few years (an incomplete 
list of references is Q— 1|. For a short state-of-the-art review see e.g. (!(]]). We therefore think 
that a model study based on a simplified nuclear force which, however, captures the essential 
features of the nucleon-nucleon interaction (long/intermediate-range attraction and short-range 
repulsion) can provide useful insights. 

In this study, we investigate the two-nucleon system in a space of momenta whose values 
are below a given cut-off A. For that, we start with a model nucleon-nucleon (NN) interaction 
given in full momentum space. The potential consists of two terms, an attractive one due to 
the exchange of a light meson with mass /i£ and a repulsive term parametrized in terms of a 
heavy meson exchange with a mass fin- The strengths parameters accompanying these terms 
are determined by fitting physical properties like e.g. the deuteron binding energy. We then 
divide the momentum space into two subspaces, spanning the values from zero to A and from A 
to oo, respectively. Consequently, the two-nucleon Hamiltonian can be regarded as a two-by- 
two matrix connecting the two momentum subspaces. By an unitary transformation it can be 
block-diagonalized decoupling the two subspaces. In this manner one can construct an effective 
Hamiltonian acting only in the low momentum subspace. The so constructed effective Hamil- 
tonian comprises the full physics for low-lying bound and scattering states with appropriate 
boundary conditions. Specifically, for the scattering states the initial momenta should belong to 
the low momentum subspace. Nevertheless, and this is an important remark, the physics of the 
high momenta is by the very construction included in the effective low momentum theory. Like 
in treatments performed in the context of chiral perturbation theory all or some of the resulting 
effective interaction can be cast into the form of a string of contact interactions of increasing 
powers in the momenta. For a consistent power counting to emerge, the coefficients accompa- 
nying these terms should be of natural size. In other words there should be a momentum scale 
(the naturalness scale A sca i e ) such that the properly normalized coefficients are of order one. 
In our exact effective theory, we can in fact precisely calculate these coefficients and check the 
naturalness property. Stated differently, the effective theory^ obtained by the exact momen- 
tum space projection plays the role of QCD and the expansion in terms of contact interactions 
for the heavy meson exchange is our model of the effective field theory. Note that throughout 
this investigation, we keep the light meson (pion) exchange explicitly. Further questions that 
can be addressed are: a) can one find a subdivision of the full momentum range such that the 
effective interaction is only weakly dependent on the precise choice of the value A, and b) is there 
relation between A and the natural momentum scale? We can also investigate the convergence 
properties of the expansion in terms of local operators since we have the exact solution to the 
problem under control. The effective interaction is naturally constructed in momentum space, 
but one can consider it in configuration space as well. Clearly, one has to expect the effective 
potential to look totally different compared to the original potential and as a consequence, the 
deuteron wave function will also change. More precisely, the projection of the original potential 



* 5 We stress that the effective theory considered here is not an effective fieid theory (EFT). This important 
distinction shouid be kept in mind. 
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into the subspace of small momenta induces non-localities in momentum space which in turn 
lead to very complicated co-ordinate space expressions. We will show some instructive examples 
of this phenomenon. 

In section II we formulate the model and the way the block-diagonalization is performed. We 



use an unitary operator proposed long time ago by Okubo Jllj] (see also ref. |12[), which leads 
to a nonlinear decoupling equation. We show that the T-matrix evaluated with the help of the 
effective interaction is exactly equal to the original T-matrix for the underlying potential. This 
is related to the fact that the component of the transformed wave function in the subspace of 
the high momenta above A is exactly zero for the low energy scattering states in which we are 



interested. This proof is different from the one given in refs.[ll] |T^]. An alternative way, which, 
however, assumes the knowledge of the the scattering states, leads to a linear equation [13]. This 
is briefly discussed in the end of the section. 

The numerical investigations are described in section III. It is shown how the nonlinear 
decoupling equation can be solved by a suitably chosen iteration. For that, one has to modify 
the original potential in order to avoid difficulties when the momentum exactly equals the 
chosen cut-off value. Physical observables, however, do not depend on this modification. We 
demonstrate (numerically) the exact agreement for the binding energy and scattering phase 
shifts using the effective and the original Hamiltonian, respectively. We then expand the heavy 
meson exchange term in a series of local contact terms with increasing dimension. This allows us 
to extract information on the naturalness property of the expansion coefficients in the effective 
interaction. As an application we calculate the expectation values of the various terms in the 
effective potential with respect to the effective deuteron wave function and low energy scattering 
states. This allows to study the convergence properties of the contact term expansion. We 
discuss the relevance of these results for a systematic effective field theory treatment of the 
two-nucleon system. We also regard the configuration space representation of the effective 
interaction. Note that some of these results were already discussed in the letter [j[q| . 

We summarize and conclude in section IV. 



2 Projection formalism for the effective theory 

In this section we develop in detail the formalism which allows to study the nucleon-nucleon 
interaction in a Hilbert space of momenta below a chosen cut-off in momentum space. The 
starting point is a given potential which involves all momentum scales and reproduces quali- 
tatively the two-nucleon phenomenology (in the S-waves). Besides being interesting in itself, 
such an effective theory in a space of low momenta can also be used to study various aspects of 
effective field theory approaches to the nucleon-nucleon interaction. For simplicity, we restrict 
ourselves to the S-waves. The formalism is, however, more general and can straightforwardly 
be extended to more complicated potentials. 

To be specific, consider a momentum space Hamiltonian for the two-nucleon system of the 
form 

H(p,p') = H (p)8(p-p') + V(p,p') , (2.1) 

where Hq stands for the kinetic energy and V for the spin-independent model force. We intro- 
duce the projection operators 

\p\ < A , (2.2) 
\p I > A , (2.3) 



V 
A 



d 3 p\p) (p\ , 
d 3 p \p) (p\ , 
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where A is a momentum cut-off which separates the low from the high momentum region. Its 
precise value will be given below. Apparently rj 2 = r/, X 2 = A, r]X = At? = and 77 + A = 1. Using 
the 77 and A projectors, the Schrodinger equation takes the form 



r/Hr] rjHX 
XHt] XHX 





(2.4) 



Obviously, the low and the high momentum components are coupled. Our aim is to derive 
an effective Hamiltonian acting on low momentum states only and which furthermore incorpo- 
rates all the physics related to the possible bound and scattering states with initial asymptotic 
momenta from the rj states. This can be accomplished by an unitary transformation U 

H — ► H' = U ] HU (2.5) 

such that 

7]H'X = XH'r] = . (2.6) 
We choose a parametrization of U given by Okubo lilt] 1 



, . (1 + AU)-V2 _ A t (1 + A4 t)-l/2 

1 A(l + A^A)" 1 / 2 (l + AAt)-V2 1 
where A has to satisfy the condition 

A = XArj . (2i 



It is then straightforward to recast the conditions eq.( |2.6| ) in a different form, 

A (H - [A, H] - AHA) rj = . (2.9) 
This is a nonlinear equation for the operator A, which takes the explicit form 

V(p, q) - J d 3 q' A(p, q')V(q', q) + J d 3 p' V(p, p')A(p>, q) 

- J d 3 q'd 3 p'A(p,q')V(q',p>)A(p',q) 

= (E q — E p ) A(p, q) . (2.10) 

Here we denoted the momenta from the r\ and A-spaces by q and p, respectively, and E q , E p 
stand for the corresponding kinetic energies. Once A and thus U have been determined, the 
effective Hamiltonian in the 77-space takes the form 

r]H'r] = 7/(1 + A^A)- 1 / 2 (H + A^H + HA + A^HA) (1 + A^A)- 1 / 2 ^ . (2.11) 



This interaction is by its very construction energy-independent and hermitean ||11|| jl2||P], 

After this block-diagonalization, the Schrodinger equation separates into two effective equa- 
tions in the respective subspaces. According to eq.(|2.5|) the connection between the eigenstates 
of H and H' is 

= C/% , (2.12) 
so that a priori the transformed problem separates into 

rjH'rj^' = Erj^' , (2.13) 
XH'X^' = EX^' . (2.14) 
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However, one of the components of has to be identically zero Q]. If both 77 and A components 
of VP' would be different from zero, one could form any linear combination thereof and would end 
up in the original space with an infinite degeneracy, which is incorrect. It remains to be clarified 
in which space the transformed scattering and bound states reside. Let us first regard scattering 
states. In which space lies if the scattering process is initiated by asymptotic momenta from 
the low momentum range spanned by 77? For that, we define the effective potential by 

V' = H'-H , (2.15) 

which is, of course, block-diagonal, since Hq and H' have this property. Note further that the 
resolvent operator of the transformed Hamiltonian H' is block-diagonal as well: 

U\z - H)~ l U = (z- H')- 1 = G'(z) =\ {Z ~ r,H ' r,) ~ 1 ° I • (2.16) 

^ (z- XH'X)- 1 J 

Obviously, z should not be in the spectrum of H. This block-diagonal form (2.16) has the 
immediate consequence that the scattering states |^ ), defined via 

\% {+) ) = limieG'(E q + ie) \q) , (2.17) 

lie in the 77 (A)-space, when the corresponding asymptotic momentum q belongs to the 77 (A)- 
space, respectively. Here, \q) is a momentum eigenstate and Hq \q) = E q \q). Immediately, the 
question arises what is the connection between the scattering states l^t 1 "^) = lim e _>o ieG{E q + 

ie) \q) and It is not obvious that eq. ( |2.12 ) holds also for these scattering states, which 

are defined through specific boundary conditions. We sketch now a proof, showing that eq. ( 2.12j ) 
is indeed valid in this case and that the relation 

= C/t|^ri+)) (2.18) 

is satisfied. The following steps appear highly plausible but do not replace a mathematically 
rigorous proof. Consider the left hand side of this equation: 

|^ (+) > = limieG'(E q + ie) \q) = limietf G(E q +ie)U\q) 

= U j liinieG(E q + ie) (^(ri + XA^ (l + A j Ay 1/2 (2.19) 
+ (A-77^A) (l + AA^y l/2 \ \q) . 
In the last line we used eq. ( |2.7| ). Let us define the operators B and C by 

B = r/Br] = (77 + A ] A) ~ 1/2 - 77 , (2.20) 

C = ACA= (A + AA^~ 1/2 - A . (2.21) 
Then eq. (|2.19j ) takes the form 

|^ (+) ) = C/t i im i e G{E q + ie) (\q) + F \q)) (2.22) 

e— »0 

with 

F = B + C + A{r] + B)- At (A + C) . (2.23) 
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Now the operator A and as a consequence B and C depend on the interaction V and are linear in 
V in lowest order. This is obvious from eq,( 2.10| ). Further it can be excluded that the state F\q ) 
is an exact continuum eigenstate of H to the energy E q , which would be required to generate a 
pole oc l/(ie) in the application of G{E q + ie). Only then the factor ie could be cancelled. As a 
consequence one obtains 

lim ieG(E q + ie) F\q) = , (2.24) 

and thus eq. ( |2.18| ) is indeed satisfied. One may also use a perturbative argument and expand 
the full resolvent operator 

limieG{E q +ie)F\q) = lim iely {Go{E q + ie)Vy\Go{E q + ie)F\q) . (2.25) 

Because of the above mentioned property of F it is clear that (p\F\q) does not contain a term 
proportional to 5(p — q) and therefore one can not generate a pole term 1 / (ze) through the 
application of Go(E q + ie) onto F\q). Consequently, each single term of the series on the right 



hand side of this equation equals to zero, which leads again to (2.24). We thus have shown 



that the scattering states initiated by momenta q from the low momentum 77-space are 

transformed into , which obey eq.( |2.13|) . The corresponding A components of ^^ + ^ are 

strictly zero. 

After considering the scattering states, we now turn our attention to the bound states. Again, 
the obvious question is: Where do the transformed bound states of H reside? If the cut-off A is 
sufficiently large then A goes to zero. This is a simple consequence of the fact that V is assumed 
to fall off sufficiently fast for high momenta. Consequently XH'X is approximately equal to 
XHX which contains only a small portion of V and thus can not support bound states at all. 



The transformed bound states have therefore to be solutions of eq.(2.13). For the physically 
reasonable choices we used in section 3 this turns out to be true. It is not known to us of which 
value of A this argument breaks down. In section 3 we shall use an iteration procedure to solve 
the nonlinear equation ( |2.10| ). It is known from IB] that the iteration converges only if the lowest 
eigenvalue of XHX is greater than the largest eigenvalue of the effective rj space Hamiltonian 
7]H'r/. If one chooses a sufficient small cut-off A, than we can expect that XHX will provide 
a bound state, which of course is not the physical one. In such a case the iteration method 



has to break down and the nonlinear equation ( 2.10 ) has to be solved in another manner. It is 



conceivable that then the true physical transformed bound state resides in the effective A-space. 
One can argue just in the same way to see the validity of the equation ( p. 12 ) for the states 



|*^ (_) ) = lim e ^ ieG(E q - ie) \q) and \V'S } ) = lim e ^ ieG'(E q - ie) \q). Therefore, the S- 
matrices in the original and transformed problem are the same: 

s& = <*f~H' (+) > = (^-wi v +) > = <* g - H i V +) > = v • ( 2 - 26 ) 

As a consequence the on-shell T-matrix element evaluated by means of the Lippman-Schwinger 
(LS) equation 

T' = V + V'GqT' , (2.27) 
yields exactly the same matrix elements as gained via the LS equation to the original problem 

T = V + VG T . (2.28) 

Note that in eq.(2.28) one integrates over the whole (infinite) momentum range whereas in 



eq.( |2.27 ) only momenta up to the cut-off A are involved. 
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An important observation is that in the subspace of momenta below the cut-off most local 
operators are non-local. For an arbitrary local operator 0(pi,p*2) = 0(pi) <5 3 (pi — P2) one 
obtains in the transformed space 



d 3 p ; U^p 1 ,p')0(p')U(p',p 2 ) 



(2.29) 



which in general contains the usual 5-function part but in addition also a strong non-local 
piece. These non-localities, which are easy to handle, are nothing but the trace of the high 
momentum components from the full space. Note, however, that the free particle Hamiltonian 
Hq and as a consequence the momentum operator p of a particle remain unchanged. Certainly, 
one could also unitarily transform the operators Hq and p. This would not yield any new 
aspects since the original and the unitarily transformed problems would be trivially identical, 
the unitary transformation only leads to change of basis. We shall not consider this trivial 
case any further. Physical observables such as the deuteron binding energy and phase shifts are 
identical in the original and transformed problems as shown before. We remark that for instance 
the average momentum in the deuteron will change in the transformed problem because the 
momentum operator is not unitarily transformed. This is not a problem since it is not direct 
observable. However, all current operators have to be unitarily transformed and this guarantees 
the equivalence of all observables. 

To end this section, we discuss an alternative way of determining the operator A, as exhibited 



in ref.[13|. That method uses the knowledge of the scattering states to the original Hamiltonian 
H. Let us now look in some more detail at this formalism. As it was already pointed out above, 
the connection between the scattering states in the original and transformed problem is given 
by 



(+) 



U 



(+) 



(2.30) 



N 77 + XAf]) (l + AU) 1/2 j] + (A - 77AU) (l + AAt 



-1/2 



(+) 



For momenta q from the 77-space this equation takes a simpler form 



-1/2 



V- (+) 



(2.31) 



since the resolvent operator of the transformed Hamiltonian H' is block-diagonal as already 
pointed out before, cf eq. ( |2.16 ). Using the trivial relation (77 + XAr])(rj + XArj) = 77 + XArj we 
conclude from the eq.( 2.3l[) for all momenta q € rj that 



A 



(+) 

?6l 



AAry 



(+) 

q G 1) 



Projecting this equation on to the state (p |, p > A, and making use of the relation 



P2 1 



$(pi -P2) + 



T(pi,P2,E, 



E, 



P2 



VI! 

E pi + ie 



(2.32) 



(2.33) 



where T(pi,p2, z) is the usual off-shell T-matrix, we end up with the following linear integral 
equation for the operator A: 



A(p,q) 



:/'(/?• '7- f d3qf A(p,q')T(q',q,E q ) 



E n — E n 



E n 



E q i + ie 



(2.34) 



->q J ^q 

Here the integration over q' goes from to A. Note this is not a usual Lippmann-Schwinger 
equation, since the position of the pole, E q , in the integration over q' is not fixed but moves 
with q. It is the second argument in A which varies, the first one is a parameter for the integral 
equation. 
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3 Applications 



This section is split into various paragraphs. First, we discuss in detail how to determine the 
operator A for a given realistic potential. We then compare the results obtained in the space 
of momenta below the chosen cut-off with the ones in the unrestricted Hilbert space. Since 
the projection formalism is exact, we can recover all physical quantities obtained from the full 
potential to an arbitrary accuracy in the smaller space with low momenta only. This extends 



the results exhibited in the letter 16]. Third, we use this model to study the expansion of 
short-range physics in terms of contact terms and draw some conclusions for the application of 
chiral effective field theories. All these calculations are naturally done in momentum space. For 
illustration, we finally show some results in coordinate space. 

3.1 Basic model and determination of the operator A 

Our starting point is a model potential which captures the essential features of the nucleon- 
nucleon (NN) interaction. We choose a momentum space NN potential with an attractive and a 
repulsive part corresponding to the exchange of a light and a heavy scalar meson, their masses 
denoted by hl and fin, respectively 

v <*'» = sk(i^£-r^)' <"> 

with t = (qi — q2) 2 - The strengths of the meson exchanges parametrized by Vl and Vu, respec- 
tively, will be determined later. This potential still contains all partial waves. The numerical 
investigations will be restricted to NN S-waves. One can work out the corresponding S-wave 
potential in closed form, 

viguq2) = ( Vh b ( - v t m ( < n+ «i;:f i ) , (3.2) 



2vrgig 2 \ \(qi- q2) 2 + V 2 H/ 

with qi 2 = \q\,2\- The nonlinear equation ( |2.10| ) simplifies correspondingly. Still, it can only be 
solved numerically. We do this by iteration starting with 

A V{p,q) 



Eq — EjJ 



(3.3) 



After four iterations we perform an average over the values of the operator A with suitably 
chosen weight factors. The choice of these weight factors is responsible for the quick convergence 
of this iteration method. A typical total number of iterations is 40 to achieve an accuracy 
of 0.0001 GeV~ 3 for the function A(p,q). The function (E q — E p ) in eq. ( [2.10D multiplying 
A(p, q) requires some caution when p and q go to A. We proceed by regularizing the original 
potential V(k',k). We multiply it with some smooth functions f(k') and f(k) which are zero 
in a narrow neighborhood of the points k' = A and k = A and one elsewhere. The precise form 
of this regularization does in fact not matter as already argued in |16|| . However, for the actual 
calculations presented here, we choose 

f(k) = 1, for KA-a and k > A + a , 

f(k) = - fl + cos MjizA+gm , for K-a<k<K-a + b , (3.4) 



2 V V b 

f(k) = - ( 1 + cos / ^(fc- A ~°) ji for A + a- 6</c<A + a, 



2 

f(k) = 0, for A-a + KKA + a-6 



S 




Figure 1: Modification of the potential due to the regularization. 



with a and b parameters of dimension [energy]. This modification of the potential V is depicted 
in fig. [l| Now A(p, q) based on that modified potential is well defined for p,q — > A. To further 
quantify the effect of this regularization, we set for simplicity 6 = 0. Then it is easy to estimate 
the difference of the regularized on-shell T-matrix, T I:eg (q,q), to the unregularized (exact) one, 
T(q,q), by 

T g {q,q) =T{q,q) - a + 0{a ) , (3.5) 

Eq — 

for all q < A — a and q > A + a. Similarly, the modification of the deuteron binding energy can 
be obtained by calculating the expectation value of (^d\V — V Ies \^D)' 

E reg = E -Aa^ D (A)A 2 k 2 dk V(A, k)V D (k>) + 0(a 2 ) , (3.6) 

Jo 

where ^o(k) is the momentum space deuteron wave function. Thus for any fixed q < A — a 
and infinitesimal e > one can choose a such that \T reg (q, q) — T(q, q)\ < e and |£' reg — E\ < e. 
Later we shall give numerical examples of the effects caused by that potential modification in 
some specific cases. We show later numerically that the effective potential V'(q,q') is affected 
only within the width a for q, q' — * A. There both V 1 and V go to zero and have no effect on 
the observables. 



We now have to fix the parameters for the potential V. We choose these as given in ref. [19|: 
V H = 7.291, V L = 3.177, fi H = 613.7 MeV and fi L = 305.9 MeV.0 That potential supports 
one bound state at E = —2.23 MeV and leads to S-wave phase shifts in the 3 Si partial wave in 
fair agreement with results from NN partial wave analysis. Thus this potential captures some 
essential features of the NN interaction. Next, we have to select a value for the cut-off A. In 



#6 For the light meson, we could have chosen the pion mass. However, since nuclear binding is largely due to 
correlated two-pion exchange, a somewhat larger value was chosen. All conclusions drawn in what follows are, 
however, invariant under the precise choice of this number. 
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Figure 2: The function A(p, q) for A = 400 MeV. 



principle, A could take any value, in particular also above the larger of the two effective meson 
masses in the potential. We shall comment on that below. We are, however, mostly interested 
in an effective theory with small momenta only and thus we start by setting A = 400 MeV. To 
be more precise, by "small" we mean a scale which is below the mass of the exchanged heavy 
particle so that one can consider the situation with a propagating light meson and the heavy 
meson integrated out and substituted by a string of local contact terms (as will be discussed in 
more detail below). 

The nonlinear equation is discretised using Gauss-Legendre quadrature points. Here we have 
choosen a = 20 keV, b = 10 keV and 100 Gauss-Legendre points. The shift in the binding 
energy due to the regularization is 0.012 KeV, which is a 0.01 permille effect. Note that it can 
be made smaller if so desired. The resulting A(p, q) is shown in fig. ||. Here we have not shown 
the function A(p, q) in the region of regularization, i.e. for p < A + a, q > A — a. In fact in this 
region the function A(p, q) goes smoothly to zero. This is guaranteed by the regularization and 
can be seen explicitly from eq. fl2.10p . 



3.2 Effective potential, scattering and bound states 

Having calculated the operator A, we are now in the position to consider observables. First, we 



need the transformed Hamiltonian H'. The determination of H' according to eq.( 2.11 ) requires 
the calculation of (r/ + A^ A)^ 1 / 2 . This is done in the following manner: as already described in 
section 2, we first define the function B(q,q') as given in eq, (|2.20| ). Consequently, this function 
has to satisfy the following nonlinear equation as follows from simple algebra: 

1 f°° 1 

B(q,q') = ~- p 2 dpA(p,q)A(p,q')-- q 2 dq B(q,q)B(q,q') (3.7) 

2 J a 2 Jo 

- J™ p 2 dp £ q 2 dqA(p, q)A(p, q) (b(q, q') + £ q' 2 dq' B(q, q')B(q', q')) . 
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The function B(q, q') denned in eq.( [2.20 ) is obviously symmetric in the arguments q, q': B(q, q') = 
B(q', q). We have solved equation ( |3,7| ) by iteration, using the same Gauss-Legendre quadrature 
points as discussed before and starting with B = — (1/2)^4^. We have found a very fast con- 
vergence of the iteration method in this case. The integrations present in eq.( p.ll ) to determine 
H' are performed by standard Gauss-Legendre quadratures and we end up with an effective 
potential V'(q',q) defined for q, q' < A. It is displayed in fig. || in comparison to the original 
underlying potential for A = 400 MeV. Note that the region of the regularization is again not 




Figure 3: Effective two-nucleon potential (solid green lines) in comparison with the original 
potential (dashed blue lines) for momenta less than 400 MeV. 

shown to keep the presentation clearer. One finds that the effective and the original potentials 
have a similar shape for momenta below the cut-off. The main effect of integrating out of high 
momentum components at the level of the potential seems to be given in this case just by an 
overall shift. 

The solution of the effective LS equation ( p. 27 ) is now very simple since the integration is 



confined to q < A. The effective bound state wave function obeys a corresponding homogeneous 
integral equation 

-^tt(g) + / q 2 dqV'(q,q)V(q) = E*(q) , (3.8) 
m N Jo 

with m^v = 938.9 MeV the nucleon mass. Using 40 quadrature points the resulting binding en- 
ergy agrees within 9 digits with the result gained from the corresponding homogeneous equation 
driven by the original potential V and defined in the whole momentum range. Furthermore, 



the S — wave phase shifts agree perfectly solving either eq.( 2.28 ) in the full momentum space 
or eq.( |2.27| ) in the space of low momenta only. This is shown in fig. ||. Note that due to 
the regularization the phase shifts go to zero for q — > A and that in the figure, the phase 
shifts are shown as a function of the kinetic energy in the lab frame and the zero occurs at 
T lab = 2A 2 /m N = 341 (85) MeV for A = 400 (200) MeV. 

One can repeat this numerical exercise choosing other cut-off values. Setting for instance 
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Figure 4: Phase shifts from the effective potential (red dashed line) and the original potential 
(green solid line) as a function of the kinetic energy in the lab frame. Left (right) panel: A =400 
(200) MeV. 



A = 2 GeV, the corresponding function A(p, q) is shown in fig. || and the effective potential 
V{q',q) is rather close to the original one, [V(q',q) — V'(q', q)]/V(0, 0) ~ 0.02. Here we have 
choosen a =200 keV and b =100 keV, which leads to the same shift in the deuteron binding 
energy as in the case with A =400 MeV. The closeness of the effective potential with the original 
one for such a large cut-off value had to be expected, since the higher momenta play only a 
perturbative role. On the other hand for rather small cut-off values of A, like e.g. 200 MeV, 
the effective two-nucleon potential is quite different from the original one. This is shown in 
fig. ^[ Again for A = 200 MeV the binding energy and the phase shifts agree perfectly using the 
effective and the original formulation, see fig. (Q). We have found solutions up to A ~ 100 MeV. 
For smaller cut-off values, the effective potential develops so much structure that one would 
have to modify the method of solving the integral equation to determine A. For the reasons 
mentioned above, we do not pursue this issue here any further. All that demonstrates that the 
nonlinear decoupling equation can be solved reliably and consequently, the effective Hamiltonian 
including all the physics can be determined to arbitrary high precision. 

To the best of our knowledge this is the first time that this projection formalism has been 
applied to such a case. In the context of the nuclear many-body problem, one often uses a 
formalism which decouples two spaces, one defining a model space of low-lying states and the 
other the rest of the -particle Hilbert space (the so-called P and Q spaces). The aim is to 
derive an effective interaction acting in the model space only. This should lead to the exact 
spectrum of the low-lying states. To do that a suitable transformation is performed to decouple 
the two spaces. It leads to a decoupling equation of exactly the form given in eq.( [2,10 ), however, 
acting now in the space of N particles. This problem is very hard to solve and approximations 
are needed. In approximate schemes that nonlinear equation is often reformulated into a linear 
form on a two-body level using the exactly known interacting two-body states (some references 



are |18| |13| ). This is also a feasible way to proceed in our context as was displayed at the end of 
section 2. Because we consider two particles only, it is an exact reformulation and the operator 
A obeys the linear equation fl2.34| ). We have solved that equation for an S-wave. For that, the 
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Figure 5: The function A(p, q) for A = 2 GeV. 



function A(p,q) was expanded into cubic splines Sk(q) 0] 

A(p,q)=Y / S k (q)A(p,q k ) (3.9) 
k 

based on a set of suitably chosen gridpoints in the interval [0, A]. This procedure leads to 
a system of linear equations for the expansion coefficients A(p, q^), which is solved by standard 
methods. Consequently, the dynamical input is the S-wave t-matrix T(q',q,E q ). The resulting 
A(p, q) agrees perfectly with the one obtained by solving the nonlinear equation. 

3.3 Implications for chiral effective field theories 

Now let us establish some contact to the theory of nuclear forces in chiral perturbation theory. In 
this approach only the nucleons and lightest mesons (i.e. pions) are typically kept in the theory as 
explicit degrees of freedom. The structure of the pion-nucleon interaction is constrained by the 
spontaneously broken chiral symmetry. To take into account the effects of heavier degrees of free- 
dom such as heavier mesons and baryon resonances, contact interactions between the nucleons 
are introduced in the corresponding Lagrangian. Such contact interactions are not constrained 
by chiral symmetry. Starting from the most general Lagrangian for pions and nucleons and 
using power-counting rules one can obtain formal expressions for the effective nucleon-nucleon 
potential (see for instance ^|{), which has to be put into the Lippmann-Schwinger equation to 
generate the corresponding S-matrix. Note that regularization and renormalization are needed 
to treat the ultraviolet divergent integrals in the Lippmann-Schwinger equation. Furthermore, 
the pionic loops in the effective NN-potential should be regularized consistently with the regular- 
ization of the Lippmann-Schwinger equation. This regularization is necessarily nonperturbative: 
the presence of a low-energy bound state signals the failure of perturbation theory. Once the 
theory is regularized (for instance with some cut-off) all parameters in the Lagrangian may be 
fixed by fitting to data, like e.g. the low energy NN partial waves and deuteron properties. It is 
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Figure 6: Effective two— nucleon potential (solid green lines) in comparison with the original 
potential (dashed blue lines) for momenta less than 200 MeV. 

commonly believed that the coupling constants scaled by some effective mass parameter A sca i e 
should have the property of "naturalness", which means they should be of the order of one. 
Only then this expansion makes sense and the power-counting is self-consistent. The value of 
the scale A sca i e is obviously closely connected to the radius of convergence of this expansion. 
Let us first consider the simpler case when the pionic degrees of freedom are integrated out. In 
that case, the low-energy NN interaction can be described entirely in terms of contact terms 
and one expects the scale A sca i e to be of the order of the pion mass m n . Therefore, all physical 
parameters which describe the NN phase shifts up to center of mass momenta of the order m^, 
such as the scattering length a and the effective range r e , are expected to scale like appropriate 
inverse powers of m n . This is, however, not the case in nuclear physics: the NN scattering length 
in the 1 S'o-channel takes an unnatural large value, a = (—23.714 ± 0.013) fm ^> l/m w . As a 
consequence, the expansion in powers of small momenta breaks down at momenta much below 
the scale m n . The physics of this phenomenon is well understood (there is a virtual bound 
state very near zero energy) and amounts to a fine tuning between different terms when the 
corresponding phase shifts are calculated |2(J0. To achieve such cancellations in the effective 
field theory calculations one has to "fine tune" the parameters. We shall see below how this 
"fine tuning" works in our model. 

The situation is more complicated in the effective theory with pions since different scales 
appear explicitely. That is why it is not clear a priori what scale enters the coefficients in the 
Lagrangian. Using a modified dimensional regularization scheme and renormalization group 
equation arguments, the authors of Q have argued that A sca i e ~ 300 MeV. On the other hand 
it was argued by the Maryland group || that no useful and systematic effective field theory 
(EFT) exists for two nucleons with a finite cut-off as a regulator. A systematic EFT is to 
be understood in the sense that the contributions of the higher-order terms in the effective 
Lagrangian to the observables at low momenta (i.e. the quantum averages of such operators) 
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are small and therefore truncation of such operators at some finite order is justified. Within our 
realistic model, we can address this question in a quantitative manner as shown below. 

Let us now apply the concept of the effective theory to our case. The original potential plays 
the role of the underlying theory of NN interactions in analogy to QCD underlying the true 
EFT of NN scattering. Integrating out the momenta above some cut-off A < \xn, we arrive via 
the unitary transformation at the effective potential V'(q,q'). In analogy to the true EFT we 
decompose V into two parts: 

V'(q',q) = V^ t (q',q) + V^ ntact (q',q) . (3.10) 

Here, Montact is a string of local contact terms of increasing dimension, which is caused by the 
heavy mass particle and the high momenta p > A. The piece Might is related to the light meson 
exchange, but again it is modified by integrating out the high momenta and in general it will 
depend on the order to which we expand M/ontact- However, since we work in a model which 
serves as an illustration, we simplify the procedure and keep Might fixed as the light meson 
exchange: Might = Might , where Might denotes the second term in eq. (|3.2|) . Specifically, 

(3.11) 

with 



(3.12) 



and the superscript '(2n)' gives the chiral dimension (i.e. the number of derivatives). Thus the 
first term on the right hand side of eq. ( |3.10| ) is the purely attractive part of the original potential 
due to the light meson exchange. In this way, we have an effective theory for NN interactions, in 
which the effects of the heavy meson exchange and the high momentum components are approxi- 
mated by the series of NN contact interactions and the light mesons are treated explicitely. Note 
that the constants Cj, C[ in eq. ( 3.12j ) correspond to renormalized quantities since the effective 
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potential V'(q',q) is regularized with the sharp cut-off A. The first question we address in our 
model is: what is the value of the scale A sca i e and its relation to the cut-off A? Of course, a 
priori these two scales are not related. For the kind of questions we will address in the following, 
we can, however, derive some lose relation between these two scales as discussed below. Since 
we know V'(q',q) numerically, we can determine the constants Cj by fitting the eq.( p.ll ) to 



V(q',q) — Might (</, <?)• This is done numerically using the standard FORTRAN subroutines for 
polynomial fits to functions of one variable and taking the corresponding polynoms typically of 
order ten to eleven. Once this is done for the functions V(q, 0) and V(q,q) all constants Cj, C[ 
in eq. ( |3.12j ) can be easily evaluated. For the potential parameters used so far and the choices 
of A = 400 MeV and A = 300 MeV, the resulting constants are given in Table 1: 
Naturalness of the coupling constants Cj, C[ means that 

— a n A scale , (3.13) 



C2n+2 

where the a n are numbers of order one. Indeed, as one can read off from table |], such a common 
scale exists, namely 

A scale = 600 MeV . (3.14) 
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C 


c 2 


c± 


ci 


c 6 


c« 


A = 400 MeV 


11.23 


-32.27 


86.73 


113.9 


-232.6 


-913.6 


A = 300 MeV 


11.15 


-32.93 


85.76 


115.6 


-265.2 


-783.7 



Table 1: The values of coupling constants Cj, C[ in [GeV^ 2 *)] for two choices of the cut-off A. 



This is a reasonable value in the sense that it is very close to the mass \in of the heavy mesons, 
which is integrated out from the theory. Stated differently, the value for A sca i e agrees with 
naive expectations. In fact, in the full theory (i.e. without the projection into the space of 
small momenta only), this observation would be trivial since the expansion of the heavy meson 
propagator gives 

a- = -5-(l--5- + -z- + ...J , 3.15 

t + ^H V (J>h Mfl > 

and thus A sca i e = fijj. In our approach, however, the large momentum components are mapped 
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Figure 7: Effective two-nucleon potential V'(q',q) (solid green lines) in comparison with the 
truncated expansion (jTllD (dashed blue lines) for A = 300 MeV. 

into the Hilbert space of small momenta and thus a priori it is not obvious that A sca i e is indeed 
given by fin- As it turns out, even for cut-off values like A = 300 MeV, there is only a small 
difference between the heavy meson mass and the ensuing natural mass scale. Another obser- 
vation is that the values of the Cj, C[ depend very weakly on the concrete choice of the cut-off 
A. Clearly, for such an expansion of the heavy mass exchange in terms of contact interactions 
to make sense, A has to be chosen below fxu- Furthermore, since we explicitely keep the light 
meson, A should not be smaller than the mass hl- If one were to select such a value for the 
cut-off, one could also consider the possibility of expanding the light meson exchange in a string 
of contact terms. We do, however, not pursue this option in here. Therefore, we conclude that 
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one should set A < A sca i e but it is not possible to find a more precise relation. In fig. [?] we 
show how well the potential V'(q',q) is reproduced by the truncated expansion eq.( 3.11| ) for 
A = 300 MeV. 

We have also calculated the two-body binding energy and the phase shifts using the form 
eq.( |3.1l| ) with the constants given in table |l[ The corresponding results are shown in fig. |8| and 
in table ||. The agreement with the exact values is good. However, one sees that terms of rather 
high order should be kept in the effective potential V^ on t a ct in order to have the binding energy 
correct within a few percent. Note, however, that the value of the binding energy is unnaturally 
small on a typical hadronic scale like the pion mass or the scale of chiral symmetry breaking. The 
phase shifts are described fairly well for kinetic energies (in the lab) up to about 100 MeV if one 
retains the first three terms in the expansion eq.( 3.11[ ). For A = 400 MeV, one can not expect 
any reasonably fast convergence for the binding energy any more. This can be traced back to 
the fact that one is close to the radius of convergence for momenta close to the cut-off. More 
specifically, with q = q' = 400 MeV, the pertinent expansion parameter is (q' 2 + q 2 )/A 
The ensuing very slow convergence is exhibited in table 0. 
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7.15 
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Table 2: The values of the binding energy calculated with Vcontact, eq. ( 3.11| ), for A = 300 MeV 
(second row) and A = 400 MeV (third row). 
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Figure 8: Phase shifts from the effective potential V' (solid line) and the truncated expansion 
(3.11) as a function of the kinetic energy in the lab frame for A =300 MeV. Left (right) panel: 
the constants Cj, C[ are fitted to the effective potential (to the NN phase shift). 



Although we have found that A sca i e ~ 600 MeV and therefore the expansion of V 1 in terms of 
contact terms seems to converge for the chosen value of the cut-off A = 300 MeV, the ultimative 
test of convergence of the expansion (3.11) would be to calculate the quantum averages of 
operators V^°\ V^ 2 \ as it was proposed in j|]. This allows to directly draw 

conclusions about the size of such operators in the effective action. In fact, the potential is 
not an observable and naturalness of its coefficients can become meaningless when one itearates 
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the potential in the LS equation since in general large momenta are involved. In our case, this 
can not happen and in the worst case, the convergence for observables will be controlled by 
the parameters A/A SC£ j e multiplied by coefficients of order one (with the exception of quantites 
which require strong fine tuning). These remarks will be substantiated by the results to be 
presented below. Because of technical difficulties, the authors of ref.fU have used for this check 
the bound-state wave function in the 1 S'o-channel, obtained from the zeroth-order potential 
and without pions (in real world there is, of course, no bound state in the 1 5o~channel). This 
allows only for a very rough estimate of the size of the operators in eq. ( 3.1 1| ) . In our model we 
can perform numerically exact calculations of these quantities using not only the bound-state 
wave function but also the scattering wave functions. Using the relation eq.( 2.33| ) one obtains 
for an arbitrary operator O 



(% (+) \0\% {+) ) = (q\0\q) +h+ II + h 



(3.16) 
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The results for quantum averages of the operators V^ n \ (n = 0,1,2,3), are shown in table |3|. 
One observes a good convergence in agreement with the naive expectation. Indeed, since the 
cut-off is choosen to be A =300 MeV and the scale A sca i ~600 MeV, one expects the expansion 
parameter to be of the order A/A sca i e ~ 1/2. Such a value agrees well with the one extracted 
from the results shown in table ||. Note further that for higher energies, the convergence is 
slower, which is also rather natural. 
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deuteron 


25.94 MeV 


-4.23 MeV 


1.07 MeV 


-0.38 MeV 


E lah = 10 MeV 


31.34 GeV~ 2 


-6.19 GeV~ 2 


1.64 GeV~ 2 


-0.58 GeV~ 2 


.Blab = 50 MeV 


11.54 GeV" 2 


-3.28 GeV- 2 


1.11 GeV~ 2 


-0.44 GeV~ 2 


E Uh = 100 MeV 


7.79 GeV- 2 


-3.09 GeV~ 2 


1.40 GeV" 2 


-0.71 GeV~ 2 



Table 3: The quantum averages of the operators V^°\ V^ 2 \ and for the bound (second 
row) and the scatterng states (third to fifth rows) for A = 300 MeV. 



In the real world to which one applies the effective field theory, one does not know the true 
effective potential V and therefore also the constants Ci,C- are unknown. They have to be 
determined by fitting such a type of effective potential to the NN data, like the low energy 
phase shifts and/or the deuteron binding energy. We can perform this exercise also in our case. 
Keeping again Might o) as a seperate piece we have determined by trial and error the constants 
Co, C2,... by solving the homogeneous and inhomogeneous LS equations and truncating the 
series at various orders. The results can be summarized as follows: For A = 300 MeV we restrict 
the fit to the phase shifts to momenta below 260 MeV. In that way, we avoid the distorsion of the 
phase shifts at the edge of its kinematically allowed region in the model space. A fit including 
only the first four terms, i.e. the constants Co, C 2 , C4 and C4 gives an excellent reproduction of 
the phase shifts as shown in the right panel of fig. ||. The corresponding constants are Co = 10.74, 
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C 2 = -23.96, C 4 = 103.7 and C' A = -288.4, all (with the exception of C' A ) within 25% of the 
exact values. The deviation of the fitted value for C' 4 from the exact one reminds us of the fact 
that such fine tuning can produce sizeable uncertainties in higher orders in such type of cut-off 
schemes and thus the interpretation of such values has to be taken with some caution. The bound 
state energy, which is not fitted, turns out to be 2.27 MeV, just 2% above the exact value. If one 
includes the sixth order terms in the fit, the results for the fourth and sixth order coefficients 
become unstable. That can be traced back to the very small contribution of these terms to the 
phase shift as long as q < 260 MeV. A stable fit can be obtained by demanding that the coefficient 
obey naturalness within a certain range, say within 20% or so. Clearly, to use polynoms of such 
high order to fit the smooth phase shifts which are already very well approximated by the first 
four terms does not make much sense. It is interesting to observe that the fit itself tends to 
limit the terms in the chiral expansion giving a good description of the phase shifts and binding 
energy. This behaviour is reminiscent of the one found in the dispersion-theoretical analysis 
of the nucleon form factors, where one gets a best fit to the hadronic spectral function with a 
limited number of vector meson poles (for details, see e.g. [21] |T~ 




Figure 9: Effective two-nucleon potential (solid green lines) in comparison with the original 
potential (dashed blue lines) for momenta less than 300 MeV in the 1 So channel. 

One can also perform a similar analysis for the two nucleon 1 Sq channel, which is expected 
to be even more troublesome than the 3 Si channel for the effective field theory approach since, 
as already stated before, the scattering lenth takes an unnatural large value. The 1 So two- 
nucleon phase shift is satisfactorily reproduced by choosing the following set of parameters 



for the potential eq.Q) (TJ: V H = 7.291, V L = 2.605, fi H = 613.7 MeV and y, L = 305.9 
MeV. All parameters with the exception of Vl remain the same as for the 3 Si channel. Vi is 
choosen somewhat smaller to make the attraction weaker. With this parameters one finds for 
the scattering length a = — 23.6fm, which is rather close to the empirical value. The reason 
for this large value is the virtual bound state of almost zero energy, which we, of course, also 
recover. We have calculated the effective potential for A =300 MeV, as shown in fig. ^. Again 
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the dominant effect of integrating out the higher momenta seems to be given by an overall shift. 
Using the same decomposition as in eq. Q3.lC ), we have again expanded ^'contact into a string of 
local contact terms of increasing dimension, cf. eq.( 3.12 ). Keeping only the terms up to fourth 
order, we have fixed the values of the constants Q by fitting eq. ( 3.1 1| ) to V'(q', 
The corresponding constants are C = 10.62 GeV" 2 , C 2 = -32.11 GeV~ 4 , C. 



- Might (<?',<?)• 

85.8 GeV" 6 



and C4 = 117 GeV -6 , rather close to the constants from table |]. (The reason for this is 
that the term in the potential corresponding to the heavy meson exchange is not modified.) 
Consequently, all conclusions about the size of the scale A sca i e from the analysis of the 3 Si 
channel are valid in this case as well. The corresponding phase shifts are shown in fig. [III. Again 
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Figure 10: Phase shifts from the effective potential V' (solid line) and the truncated expansion 
eq.( |3.11| ) as a function of the kinetic energy in the lab frame for 1 S'o channel for A =300 MeV. 

as in the case of the 3 Si channel, one sees that the phase shift is described well up to about 
100 MeV if the terms up-to-and-including fourth order in the expansion e q.(|3.11| ) are kept 
explicitly. To further study the convergence of the contact term expansion eq. (|3.11|) in the 1 So 
channel, we have again calculated the expectation values of operators V(°\ and using 
the scattering states for different energies. As exhibited in table || the values of the quantum 
averages are slightly different from the ones for the 1 S'o channel and the expansion eq.( |3.11| ) is 
again convergent. However, one observes a much slower convergence for the scattering length, 
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E Uh = 10 MeV 


37.42 GeV" 2 


-4.68 GeV" 2 


0.90 GeV" 2 


.Blab = 50 MeV 


11.72 GeV" 2 


-2.56 GeV" 2 


0.68 GeV- 2 


Blab = 100 MeV 


8.82 GeV" 2 


-2.98 GeV" 2 


1.18 GeV" 2 



Table 4: Quantum averages of the operators V^°\ and for the scattering states with 
A = 300 MeV in the l S channel. 

again because of its unnatural large value. Keeping V^°\ + y( 2 ) and + + in the 
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expansion for V^ on t ac t one finds a = — 9.0fm, a = — 54.5fm and a = — 21.2fm, respectively. This 
is analogous to the situation with the binding energy in the 3 5i channel, compare table | ( some 
recent references dealing with the particular features of this partial wave in the EFT approach 
are [0 , @] ,[§[]) • Certainly, one can achieve a considerably faster convergence with respect to the 
scattering length by fine tuning the constants Cj (i.e. by fixing them from a direct fit of the 
effective potential to the phase shift), as it was illustrated for the 3 Si channel. 

3.4 Coordinate space representation 

As discussed before, applying the method of unitary transformation turns local operators in 
momentum space into non-local ones. Similar effects happen in the coordinate space represen- 
tation. To make this transparent, we consider here the configuration space representation of the 
effective potential V . Denote by Vi(p',p) the effective momentum space potential for angular 
momentum I. The corresponding r-space expression is obtained from 

2 f 

Vi(x',x) = - p 2 dpp' 2 dp' ji(px)Vi(p',p)ji(p'x') , (3.18) 
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with ji(y) the conventional I spherical Bessel function. The S-wave potential (I = 0) for 



A = 400 MeV is shown in the left panel of fig. 11. It is, of course, symmetric under the 



interchange of x and x', but looks very different from a local potential. 




Figure 11: Coordinate space representation of the effective S-wave potential Vo(x,x') . Left 
panel: A = 400 MeV and < x,x' < 6fm. Right panel: A = 5.5 GeV and 1 < x,x' < 2fm. 

Of course, if one increases the value of the the cut-off A, a peak along the diagonal x = x' 
resembling the delta function should develop and the superposition of the two Yukawa potentials 
related to the light and heavy meson exchanges, respectively, should appear along the diagonal. 



This is indeed the case as demonstrated in the right panel of fig. 11 for A = 5.5 GeV. Note also 
that in this case, where U ~ 1, the range of the potential is essentially given by the inverse of 
the light meson mass. One can also construct the momentum and coordinate representations of 
the deuteron wave function. For a momentum space picture, we refer to ref. |l 
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4 Summary and conclusions 



In this paper, we have shown how to construct an effective low energy theory for nucleons based 
on the method of unitary transformations starting from a realistic two-nucleon potential (in 
momentum space). This unitary transformation can be parametrized by an operator A, which 
obeys a nonlinear integral equation. This equation can be solved numerically and any observable 
can then be calculated in the space of small momenta only. To the best of our knowledge such an 
exact momentum space projection has never been done before. While the method is interesting 
per se, we have also made contact to chiral perturbation theory (CHPT) approaches to the 
two-nucleon system by studying a series of questions, which can be addressed unambigouosly 
within the framework of our exact low momentum theory. Clearly, this should not be considered 
a substitute for a realistic CHPT calculation but might be used as a guide. The salient results 
of this investigation can be summarized as follows: 

1) We have demonstrated that the theory projected onto the subspace of momenta below a 
given momentum space cut-off A leads to exactly the same S-matrix as the original theory 
in the full (unrestricted) momentum space provided appropriate boundary conditions for 
the scattering states are chosen. In particular, the components of the transformed scat- 
tering states with initial momenta below the cut-off A in the subspace of momenta above 
the cut-off A are strictly zero. It is important to stress that the exact projection leads to 
non-localities in momentum space. 

2) Starting from an S-wave NN potential with an attractive light (hl ~ 300 MeV) and 
repulsive heavy meson exchange (fj,jj — 600 MeV), we have numerically solved without 
any approximation the nonlinear equation for the operator A and demonstrated that the 
bound and scattering state spectrum of the effective and the full theory agree exactly up 
to the cut-off A. In particular, we have exactly one bound state with a binding energy of 
2.23 MeV. These results are independent of the value of the cut-off, which was varied from 
200 MeV to 5.5 GeV. We have argued that the most natural choice is A about 300 MeV. 
The effective potential can differ substantially from the original one (for values of A on 
the small side of the range mentioned before). 

3) We have expanded the heavy meson exchange term in a string of local operators with 
increasing dimension but kept the light meson exchange explicitely. The corresponding 
coupling constants accompanying these local operators, which are monomials of even power 
in the momenta, can be determined precisely from the exact solution. We have shown that 
they are of "natural" size, i.e. of order one, with respect to the mass scale A sca i e = 600 MeV. 
We have also discussed the relation of this scale to the mass of the heavy meson, which is 
integrated out, and the convergence properties of such type of expansion. In particular, 
to recover the binding energy within a few percent, one has to retain terms of rather 
high order in this expansion, cf. eq.( |3.11 ). This is to be expected due to the unnatural 



smallness of this energy on any hadronic mass scale. The 3 Si scattering phase shift can 
be well reproduced up to kinetic energies T\ a \> ~ 120 MeV with the first three terms in the 
contact term expansion. 

4) Based on the expanded heavy meson exchange term, we have also determined the constants 
Ci directly from a fit to the phase shifts. This is equivalent to the procedure performed 
in an effective field theory approach. We could show that as long as one does not include 
polynoms of order six (or higher), the resulting values of these constants are close to 
their exact ones. Furthermore, the binding energy is reproduced within 2%. Including 
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dimension six terms, the fits become unstable. This can be traced back to the fact that 
the contribution of such terms to the phase shifts are very small (at low and moderate 
energies) and thus can not really be pinned down. 

5) We have also studied the quantum averages of the expanded potential in the bound and 
scattering states. For A = 300 MeV, the expansion parameter is of the order of 1/2 and 
we find fast convergence for the bound and the low-lying scattering states. As expected, 
for scattering states with higher energy, the convergence becomes slower. 

6) To study the 1 Sq channel, we had to slightly readjust the parameters of the model potential. 
The phase shift can be well reproduced with the terms up-to-and-including fourth order 
in the contact term expansion of the heavy meson exchange. For the scattering states, 
the quantum averages of the expanded potential show convergence properties similar to 
the 3 Si case. There is no bound state in the 1 Sq, but a virtual one just above threshold. 
Therefore, the pertinent scattering length is unnaturally large and it shows a similar slow 
convergence as does the binding energy in the 3 Si channel. 

7) In the model space of small momenta only, one can also study the non-localities in the co- 
ordinate space representation. We have shown that for typical cut-off values, the effective 
potential V(x, x') is highly non-local and looks very different from the original one. For 
very large values of the cut-off, one recovers the original local potential. 

We hope that this study might be useful for derivation of NN-forces based on chiral Lagrangians 
in the low-momentum regime. It should also provide new insights into a consistent and conver- 
gent treatment of relativistic effects in few- and many-nucleon systems. 
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